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We show that a “pile up” effect occurring for a train of non-interacting identical particles incident 
on the same side of a resonance scatterer leads to significant interference effects, different from 
those observed in Hong-Ou-Mandel experiments. These include characteristic changes in the overall 
transmission rate and full counting statistics, as well as ’’bunching” and ”anti-bunching” behavior in 
the all-particles transmission channel. With several resonances involved, pseudo-resonant driving of 
the two-level system in the barrier, may also result in sharp enhancement of scattering probabilities 
for certain values of temporal delay between the particles. 
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I. INTRODUCTION 

Quantum statistical effects accompanying scattering of 
non-interacting identical particles, are among some of the 
most intriguing predictions of quantum mechanics. Their 
studies, both theoretical and experimental, now consti¬ 
tute an extensive research field m-m- If two such parti¬ 
cles, prepared in wave packet states meet head on in free 
space, they would eventually “pass through each other”, 
just like their distinguishable counterparts. The situa¬ 
tion is different if such particles coincide inside a scat¬ 
terer, with a possibility of two (or more) distinct scatter¬ 
ing outcomes for each particle. In the celebrated Hong- 
Ou-Mandel (HOM) experiment [5], the particles enter¬ 
ing a scatterer from opposite sides, are seen to leave the 
barrier predominantly from the same side (bosons), or 
from the opposite sides (fermions). The HOM effect has 
found important practical applications in quality testing 
of single-photon sources [4], entanglement detection [5], 
entanglement swapping [B], and quantum metrology [7]. 
Its generalizations, to name a few, include observation of 
multiple photon bunching effects [SHini scattering of pho¬ 
tons by multiport beam splitters and using them as in¬ 
terferometers for identical particle in spatially separated 
modes [IIHS], Efforts to extend the HOM interference 
experiments to bosononic or fermionized cold atoms m 
are currently underway [HUS]. 

Perhaps surprisingly, little studied to date remains the 
case complementary to that of Hong, Ou and Mandel, in 
which the particles enter the barrier from the same side 
and “meet” there owing to a “pile up effect”, if the bar¬ 
rier is capable of detaining the first particle long enough 
for the following particle(s) to catch up with it. The sta¬ 
tistical effects are, in this case, quite different from those 
predicted for the HOM interference, and are most pro¬ 
nounced in resonance tunneling, where a particle spends 
in the barrier roughly the lifetime of the corresponding 
metastable state, which can be large for sharp resonances. 

In this paper, we give the general theory of the ef¬ 
fect and demonstrate how the said “piling up” alters the 
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FIG. 1. (Color online) Schematics diagram showing a ’’train” 
of N wave packets incident on a double barrier which supports 
two resonance levels. 


transmission rate for initially correlated many-particle 
states. There are complex interference effects in the scat¬ 
tering statistics of the incident particles, whose wave 
packet modes do not overlap prior to their arrival at 
the scatterer. A preliminary analysis of the case of two 
fermionized atoms can be found in |24j . In |25j it was 
demonstrated that interference effects of the similar kind 
would arise whenever a particle simultaneously populates 
several wave packet modes. For brevity we will refer as 
particles (fermions or bosons) to both cold atoms and 
photons, equally amenable to our analysis. 

The rest of the paper is organised as follows: in Section 
H we discuss the correlation between initial particles. In 
Section HI we analyse the correlations acquired in scat¬ 
tering. In Section IV we introduce a generatng function 
for the scattering statistics. In Section V we show how 
quantum statistical effects vanish for particles well sepa¬ 
rated initially. Sections VI and VH discuss the two- and 
iV-particle cases, respectively. In Section VHI we aplly 
our general approach to resonance transmission across a 
scatterer supporting one or more metastable states. Sec¬ 
tion IX contains our conclusions. 








III. CORRELATIONS BETWEEN THE 
SCATTERED PARTICLES 
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II. INITIAL CORRELATIONS BETWEEN THE 
PARTICLES 

Consider, in one dimension, a source sending N identi¬ 
cal noninteracting particles of mass ^ |26j in wave packet 
states, tpn{xn), as is illustrated in Fig.l, 

'ipnixn, t) = (27r)“^/^ J An{p) exp[ipXn “ iE{p){t + tn)]dp , 

E{p)=py2p (1) 

towards a finite-width potential barrier at times 0 = < 

^2 < < t]^. 

The operators and a„, creating and annihilating an 
incident particle particle in the state V'n, are given by 

o+ = (27r)“^/^ j A{p) exp[-iE{t + t„)]a+(p)dp, (2) 

where the dagger denotes Hermitian conjugation, and 
the plane wave creation and annihilation operators 
a+(p) and a{p) obey the usual commutation relations, 
[a{p),a+{p')]^ = a{p)a+{p') =Fa+(p')a(p) = 27r(5(p-p'): 
with the upper and lower signs corresponding to bosons 
and fermions, respectively (cf. also Ref. |17]). Their 
(anti) commutators coincide with the overlaps between 
the wave packets Q, 

[am,a+]^ = J A*^{p)Anip)exp[iE{p){tm-tn)]dp 

= Imn (3) 

where Inn = 1 and Imn = Inm Ttie symmetry of the inci¬ 
dent state has no effect on the initial probability density 
provided Inm = dnm, e.g., for the delays between emis¬ 
sions, \tm — tn\, large enough for the rapid oscillations 
of the exponential in ([^ to destroy the integral for all 
m ^ n. We will refer to such particles as initially uncor¬ 
related. Alternatively, the particles may be prepared in a 
correlated initial state, and below we will consider both 
these cases. It is readily seen that spreading of freely 
moving wave packets doesn’t alter the commutation re¬ 
lations ([^. Thus, the symmetrized or anti-symmetrized 
wave function describing N incident particles is given by 

N 

\^Ut)) = n “n W|0)> (4) 

n—1 

where |0) is the vacuum state, and K is the normalization 
constant. By Wick’s theorem, we have (the upper and 
lower signs are for bosons and fermions, respectively) 

N 

K=J2 (±1)^("('^» n 4-. = [Imn] (5) 

(t{N) "i —1 

where cr{N) is a permutation of the indices (0,1,..,A^) 
and p{a) is its parity [SH] . 


At large times, after all particles have left the barrier 
area, each wave packet ends up split into the transmitted 
(t) and reflected (r) parts. Thus, as t —)■ oo, the wave 
function has the form 

N 

l'kout(i)) = K-^/^l[[tt{t) + r+{t)]\0) . (6) 

n—1 


where the corresponding creation and annihilation oper¬ 
ators are given by 

= J T{p)An{p)exp[-iE{t-\-tn)]a~^{p)dp , 

r+= y R{p)Anip)exp[-iE{t-\-tn)]a~^{-p)dp , 

tn = {t+) + , r„ = (r+)t, (7) 

and T(p) and R{p) are the barrier transmission and re¬ 
flection amplitudes for a particle with a momentum p. 
Since |T(p)p -|- |i?(p)P = 1, as t —>■ oo we also have 


Trr 


— [I'm J ] =F 

= [ \Tip)\'^^mip)^nip)exp[iE{p){tm-tn)]dp , 


Rmn — [xm ; ] =p — Imn Emn j 


( 8 ) 


while all remaining commutators vanish. In Eqs. ([^ 
Tmn = is a Hermitian matrix of the overlaps be¬ 
tween the transmitted parts of the wave packets, and its 
diagonal elements coincide with the probabilities 
for the n-th particle to be transmitted on its own, 


Tnn = J \Tip)\'^\An{p)\^dp = Wn- (9) 


We note that even initially uncorrelated particles may be¬ 
come correlated as a result of scattering. This would hap¬ 
pen, for example, if each transmitted one-particle state 
is significantly broadened in the coordinate space (nar¬ 
rowed in the momentum space), so that the integrals in 
Eq. (|^ do not vanish, even if the integrals in Eq. <§ did. 


IV. THE GENERATING FUNCTION AND 
FULL COUNTING STATISTICS 

For N identical particles, there are V -|- 1 outcomes, 
with n = 0,1,..., V particles crossing in the barrier 
whose probabilities, IT(n, A), we will study next. It is 
convenient to construct a generating function G(a), 

G^{a)= lim (Tout(i)l^(^)«)) > (10) 

t —¥00 

where |^'(t,a)) = IlLiK™ (1)+’’« (0]|0), and AT 

is defined by Eq. ([^ . By Wick’s theorem, we have 

G^ia) = S^[AiJ/S^[Imn] 


( 11 ) 


3 


where the matrix A is given by 


is given by 


= Imn + (« “ l)Tmn, n,m=l, 2, N. (12) 
For the mean number of transmissions, ut = 

J 2 n =0 

nTiN,ti,t2, ...,tAr) = daG^{a)\a=i = (13) 

N 

1=1 

where Imn is the matrix obtained from Imn by replacing 
the elements of the _)-th row, Iji,I jn with Tji,..., T^jv- 
The full counting statistics of an A-particle process are 
evaluated by noting that 

lF±(n,A,ti,t2,...,tjv) = (14) 

nl 

h<h<--<jrL 

where is the matrix obtained from Rmn by 

replacing the elements of the rows ji,j 2 , with the 

corresponding rows of the matrix Next we briefly 

discuss what would happen if the particles were not iden¬ 
tical. 


V. THE DISTINGUISHABLE PARTICLES 
LIMIT 


N 


N) = E] 'Wi = wN. 


(16) 


i=l 


It is readily seen that the DP limit ( |17[ ) is reached if iden¬ 
tical particles arrive at the scatterer after long intervals, 
|ti—td —)■ oo. Indeed, in this limit all operators in Eqs.([^ 
and commute, both Imn and Tmn are diagonal, and 


for Wi = Wj = w, Eqs. (14) yield 


lE±(n,7V,ti,t2,...,%) ^lE^^(n,A). (17) 


We are, however, more interested in the case where quan¬ 
tum statistical effects do lead to measurable changes in 
in the channel probabilities W^{n,N,ti,t 2 , and 

will consider it next. 


VI. THE TWO-PARTICLE CASE (N=2) 


In the simplest case of just two particles, N 
Eqs. (13) and (14) yield 


W=^( 2 , 2 ) 

W±(l, 2 ) 


riT 


WiW2 ± |ri2p 

l±|/l 2 P ’ 

wi(l — W 2 ) + ^ 2(1 — wi) ± 2Re(Ti2i?i2) 

l±|/l 2 P ’ 

(18) 

wi + W2± 2Re[Ti2/i2] 
l±|/l 2 P ’ 


The appearance of correlations between initially uncor¬ 
related particles can be explained in the following way. 
The particles are well separated initially, and if they leave 
the scatterer quickly enough, each scattering event occurs 
independently. If, on the other hand, the scatterer de¬ 
tains each particle for a significant period of time, Bose or 
Fermi statistical effects become important while several 
(or all) particles are still inside. This, in turn, may alter 
measurable probabilities for various outcomes, which is 
the effect we seek to describe here. This is not possible 
if the particles are distinguishable, and do not interact 
with each other. Such particles cannot ’’meet” in the 
scatterer (or anywhere else), are always transmitted in¬ 
dependently, and it does not matter if they arrive all at 
the same time, or their arrivals are separated by long 
time intervals. 

Assume, for simplicity, that the individual tunnelling 
probabilities in Eq.([^ are equal for all particles, Wi = 
Wj = w. Then the probability for n out of N distin¬ 
guishable particles (DP) to be transmitted is given by 
the binomial distribution, 

W°^{n,N) =C^w^{l-w)^-^, (15) 

where is the binomial coefficient. The mean num¬ 
ber of transmissions, also independent of the choice of 


which coincides with the results of [23] if the particles 
have the same momentum distribution, Ai{p) = A 2 {p). 
The last of the Eqs.(18) shows that if one sends to the 
scatterer initially correlated pairs of identical particles, 
(/12 7 ^ 0 ), the mean number of transmissions per pair 
may be different from the result obtained for two distin¬ 
guishable particles in the same wave packet states. 

If the pairs are not initially correlated, we always have 
Ut 7 ^ n^Px, but the bosons (fermions) are more (less) 
likely to exit the scatterer on the same side. This be¬ 
haviour will be observed if the two particles, initially well 
separated from each other, meet in the scatterer, so that 
Ti 2 7^ 0 . 


VH. THE A-PARTICLE CASE 


Equations (13) and (14) show that the results of the 


previous Section also hold for an arbitrary number of 
particles, N > 2. The mean number of transmissions 
may be affected by the symmetry of the initial state, 
Ut 7 ^ ^f=ijwj, if, and only if the particles are corre¬ 
lated initially, Imn 7 ^ ^mn- 

For initially uncorrelated particles, Imn = ^mn, the 
symmetry changes the probabilities W{n,N), but not 
n, provided Tmn 7 ^ WnSmn- In this case, “bunching” 
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and “anti-bunching” types of behavior can be observed 
for bosons and fermions in the probability for all N 
particles to be transmitted, W{N,N). Since the matrix 
Tmn is positive definite, the Hadamard inequality for 
determinants m, and its analog for permanents |32j 
ensure that W±{n, (cf. also Ref. [53]). 

Thus, N bosons (fermions) are more (less) likely to 
be transmitted all together than DPs in the same 
one-particle states. Note that this argument cannot 
be extended to the probabilities W{n < N,N), or for 
initially correlated initial states, Imn mn • 


VIII. RESONANCE TRANSMISSION 


0 0,5 1 1,5 2 
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A system likely to show these effects is a resonance 
barrier, where, due to the long delay in traversing it, 
even the particles well separated initially have a chance 
to “pile up” in the scatterer. Transmission coefficient of 
such a barrier can be written as a sum of narrow Breit- 
Wigner peaks. 


\T{p)? 


E 


r 


2 

I 


(pV2/r-if[)2+rr 


(19) 


and even for initially uncorrelated particles, the shape of 
|r(p)P|A(p)p may be narrow enough to ensure that 
is not diagonal even if Imn = ^mn- 

Figure 2 shows the mean number of transmissions for N 
particles emitted after equal intervals, \tm — tn\ = T, in 
identical Gaussian states of a coordinate width a and a 
mean momentum pq, 


FIG. 2. (Color online) Mean number of transmissions for 
N particles in identical Gaussian wave packets vs. time 
T between the emissions: (a) through a single resonance 
level El/Eo = 0.41, Fi/Eq = 0.0087, poa = 3.77, and 
Eo = Pof2fi ; (b) through two resonance levels Ei/Eo = 0.95, 
Fi/Eo = 0.038, El/Eo = 3.82, Fz/Eo = 0.28 and poo = 6.04. 
Horizontal dashed lines correspond to the limit of distinguish¬ 
able particles. The momentum distribution |A(p)p and the 
transmission coefficient |T(p)|^ are shown in the insets. 
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A{p) = An(p) = (cr^/27r)^/^exp[(p-po)^cr^/4]. (20) 

The scatterer supports two resonant metastable states 
with the energies EI 2 and widths ri _2 of which one, or 
both can be accessed by the incident particle, as shown 
in the insets in Figs. 2 a and b. With only one level 
involved, the mean number of transmissions nriT) for 
bosons raises to a maximum value for some correlated ini¬ 
tial state (illustrated in Fig. 3 together with its fermionic 
counterpart), and then returns to the DP limit for ini¬ 
tially uncorrelated particles (see Fig.2a). For fermions, 
the Pauli principle mostly reduces np to levels below the 
DP level, which for the maximally correlated states, ob¬ 
tained as T —>• 0 [3D], is considerably reduced. 

With two metastable states involved, interference be¬ 
tween resonances reverses the effect: for 0.1 < pg/2/rT < 
0.25, nr is suppressed for bosons, and enhanced for 
fermions (see Fig.2b). 

The scattering probabilities W (n, N, T) for the single res¬ 
onance case, plotted in Fig.4 for V = 4, show smooth 
deviations from the DP values in Eq. 0 before reach¬ 
ing those values as the times between arrivals tends to 
infinity. We note that the probabilities W{N,N) never 
fall below (exceed) their DP level for bosons (fermions), 
as discussed in the previous Section. We note also that 


FIG. 3. (Golor online).One-particle density p{x) (normal¬ 
ized to 2) for initially correlated 2-particle Gaussian state, 
Ai(p) = A 2 {p) ~ exp[-(p-po)cr^/2], Poo = 6, poT/po = 1.5, 
Pot/2p = 4.5, and T = t 2 — ti. Also shown by the dashed 
line is p{x) for two distinguishable particles in the same one- 
particle states. 


the increase or decrease in np results from a similar in¬ 
crease or decrease in the probability of the one-particle 
transmission channel, W{1,N). 

With two resonances accessible to the particles, the 
picture is more interesting. For bosons, W (n, N, T) ex¬ 
hibit maxima, whenever the time between the emissions 
coincides with a multiple of the difference of resonant 
energies, T « = 2Trk/{E2 — E^), k = 1,2,... The 

peaks are most pronounced for the (N,N) channel (cf. 
Fig. 5a) and, as shown in Fig.6a, become sharper as 
N increases. This is another consequence of the sym- 
metrization of the initial state which, with each particle 
distributed between the wave packets in Fig.l, appears 
to produce quasi-periodic excitation of the metastable 
two-level system supported by the barrier. With the 
number of particles increasing, the excitation looks more 
periodic, and the “resonance” condition T « T/j needs 
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FIG. 4. (Color online) a)Probabilities for n=l,2,3,4 bosons 
to be transmitted for = 4 (single resonance); b) same as 
a), but for fermions. The parameters are as in Fig. 2a. Hori¬ 
zontal dashed lines indicate the the corresponding values for 
distinguishable particles given by Eq.(17l Incident particles 
may be considered uncorrelated for EqT ^ 0.4. 
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FIG. 5. (Color online) a)Probabilities for n=l,2,3,4 bosons 
to be transmitted for N = 4 (two resonances); b) same as a), 
but for fermions. The parameters are as in Fig. 2b. Incident 
particles may be considered uncorrelated for EqT > 0.4. 


to be satisfied with ever greater accuracy. Note that a 
similar (yet not identical) interference effects have been 
predicted for scattering trains of wave packet modes 
representing a single particle (for details see [25]). 

For fermions, probing two resonance states, the peaks at 
T = Tk are replaced by dips, which appear, for example 
in the probability W{2,4) shown in Fig. 5b. In contrast 
to the bosonic case, these dips are never seen in the 
{N, N) channel, where W{N, N, T) undergoes sinusoidal 
oscillations, no matter how large is the numbers of 
particles N (see Fig. 6b). Experimental observation 
of effects of the Pauli principle on resonance tunneling 
would be possible for cold atoms in the Tonks-Girardeau 
regime, injected into quasi-one-dimensional trap with 


FIG. 6. (Golor online) a)Probabilities for all bosons to be 
transmitted for N = 2,4,6 (two resonances); b) same as a), 
but for fermions. The parameters are as in Fig. 2b. 


laser-induced barriers [Mj. An optical realization of the 
bosonic experiment would consist in sending identically 
polarized photons toward a Fabry-Perot interferometer, 
or injecting them in a waveguide with narrowing, 
imitating a one-dimensional barrier. If required, a 
correlated initial state can be produced by scattering 
several uncorrelated particles off a long-lived resonance, 
and selecting the outcome in one of the n-particle 
transmission channels. 


IX. CONCLUSIONS AND DISCUSSION 


In summary, Bose or Fermi statistics can significantly 
change the scattering outcomes for a train of non¬ 
interacting identical particles impinged on the same side 
of a scatterer. Physically, the effect requires simulta¬ 
neous presence of several particles inside the scatterer. 
This may occur if the initial state is already correlated, 
or if the particles, well separated initially, ’’pile up” in¬ 
side, as a result of a scattering delay. In the latter case, 
the statistics of the process are affected in such a way 
that the mean number of transmissions per train, ny, re¬ 
mains unaffected. For a correlated initial state, ut may 
be larger or smaller than that for the train composed of 
distinguishable non-interacting particles. 

Mathematically, this is an interference phenomenon aris¬ 
ing from the presence of additional terms in the (anti-) 
symmetrized wave function, which disappears if the par¬ 
ticles can be distinguished. Its analysis is extremely sim¬ 
ple, owing to the commutation of the evolution and sym- 
metrisation operators: it is sufficient to first solve the 
corresponding one-particle problems, and then evaluate 
the asymptotic exchange integrals as t —>■ oo. Interference 
plays the most fundamental role in quantum mechanics, 
and we think it unlikely that a more detailed or more 
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’’physical” explanation of the effect can be provided. 

To conclude, scattering of trains of identical particles of¬ 
fers a variety of interference effects, very different from 
those observed in the HOM experiments, some of which 
were discussed in detail in Sect. VIII. Observation of 
such effects is within the capability of modern experi¬ 
mental techniques. 
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